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\^ , Abstract. In analogy with the standard derivation of the Schwarzschild solution, 

CS| ' we find all static, cylindrically symmetric solutions of the Einstein field equations for 

vacuum. These include not only the well known cone solution, which is locally flat, 
O ! but others in which the metric coefficients are powers of the radial coordinate and the 

space-time is curved. These solutions appear in the literature, but in different forms. 



^ . corresponding to different definitions of the radial coordinate. Because all the vacuum 

solutions are singular on the axis, we attempt to match them to "interior" solutions 



> 



X 



with nonvanishing energy density and pressure. In addition to the well known "cosmic 
string" solution joining on to the cone, we find some numerical solutions that join on 
QQ ' to the other exterior solutions. 

\o ' 

PACS numbers: 04.20.-q 
1. Introduction 

Static solutions of Einstein's equations with spherical symmetry (the exterior and 
interior Schwarzschild solutions) are staples of courses in general relativity. Solutions 
with cylindrical symmetry (combining translation along and rotation around an axis) 
are much less familiar. In this paper we construct and study all the vacuum solutions 
of this sort, and we set up the differential equations for nonvacuum solutions and 
find some solutions numerically. The article is intended for students at the advanced 
undergraduate or beginning graduate level; instructors of introductory general relativity 
courses may also wish to use portions of it as a basis for exercises appropriate for 
their classes. The paper is based on the undergraduate research thesis of Cynthia 
Trendafilova [T]. 

In Section 2 we write the most general expression for a space-time metric with 
static and cylindrical symmetry and solve the Einstein vacuum field equations for the 
components of the metric tensor. Equivalent solutions may look different because 
of different definitions of the radial coordinate in a cylindrical system. At least 
three different conventions ( "radial gauges" ) have been frequently used in the previous 
literature (see Section 6), none of which is the one that seemed most natural to us. 
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Namely, we adopt the convention that the angular term in the metric is d0^. Then 
we carefully remove all redundant solutions corresponding to the freedom to rescale the 
coordinates. The principal result is the general solution ( fT8|) . in which all the metric 
components are powers of r with exponents related by the constraint (fT6l) . 

In Sections 3 and 4 we explore the geometrical natures of the various solutions 
obtained. They fall into two classes, the more plausible of which shows the circumference 
27rr increasing with radial arc length (another natural radial coordinate). Also, the 
transformations into the other radial gauges are worked out; in all cases, the metric 
functions are powers of whatever radial coordinate is adopted. Our investigation of 
the relations among these various coordinate conventions is more systematic than any 
previous one known to us. 

The space-times ( JTSl) . which are curved in general, form a two-parameter family; it 
includes a one-parameter subfamily of spaces that are locally flat. When maximally 
continued down to a central singularity, such a space-time is a cone formed from 
Minkowski space by removing a wedge, characterized by a deficit angle. Although the 
other solutions were discovered in the early days of general relativity I2l[3], they are 
much less well known than the cones. 

Like the exterior Schwarzschild solution (when it is not treated as a black hole), one 
expects vacuum cylindrical solutions to be physically relevant only over some subinterval 
of the r axis. In Section 5 we search for nonsingular "interior" solutions, with nonvacuum 
sources, that can be joined on to the vacuum solutions at some radius. The easiest to 
find have a locally flat cone solution on the outside and a constant energy density p 
inside, with pressure Pz = —p along the axis and vanishing pressures Pr = p^j, = in the 
perpendicular plane; these are the well known "cosmic string" solutions jlllH]. Although 
natural from the point of view of gauge theory (see, e.g., |6j), such an equation of 
state would be surprising for normal matter. We write down the generic Einstein and 
conservation equations, isolate the independent dynamical variables and constraint, and 
construct some numerical solutions with isotropic pressures. These solutions match onto 
the exterior nonflat vacuum solutions. The results are examples of space-times that 
also deserve to be called "cosmic strings" but have a more general equation of state 
for the matter inside the string. We do not claim that these solutions necessarily are 
astrophysically plausible. In particular, we have not investigated their stability under 
perturbations that violate the assumed symmetries. 

In Section 6 we present some historical background on the solutions and gauge 
choices presented in this paper, while Section 7 offers some concluding remarks. 

2. Vacuum solution of Einstein equations 

In order to solve the Einstein equations for a cylindrically symmetric spacetime, we must 
first determine the line element with which we would like to work. Let t denote a time 
coordinate. For a fixed time t, a cylindrically symmetric spacetime can be described as 
follows. There is a central axis of symmetry, with z denoting the coordinate along this 
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axis. Then denotes the coordinate which measures an angle around this axis, and r is 
the radial coordinate which increases as one moves away from the axis. The line element 
ds^ is then given by ds^ = ga^dx'^dx^ , where ga0 is the metric tensor describing the 
spacetime. Since we want static, cylindrical symmetry, we require that y—ds^, the clock 
time measured by a test particle in moving from {t, r, 0, z) to {t, r, 0, z) + {dt, dr, d0, dz), 
is independent of t, 0, and z. It must also be invariant under reversals in the directions 
or signs of dt, d0, and dz. Hence the general static, cylindrically symmetric metric is of 
the form 

ds^ = -e^^dt^ + e^^dr^ + e^^d^P^ + e^^dz^, (1) 

where $, A, Q, and are functions of r only. By writing our unknown functions in 
the form of exponentials, we guarantee that our coefficients will be positive as we would 
like them to be, and also mirror the standard textbook treatment of the spherically 
symmetric metric. 

In analogy to the standard treatment of spherical symmetry [71 p 256], we define 
r so that the coefficient of d0^ is equal to r^. We shall call this convention "tangential 
gauge." In Sections 3 and 4 we discuss three alternative conventions (out of infinitely 
many possible), and also the question of whether any generality is lost by this convention. 
Thus by setting e^^ = the metric can be written 

ds' = -e^-f + e^^dr^ + r^cj)'' + e'^'dz^, (2) 

where $, A, and \l/ are the unknown functions of r for which we would like to solve. 
The form in which we have written the metric does not restrict the range of to be 
from to 27r; instead it runs from to some angle 0^,. As we shall show later, can 
be forced to fill an angle of 2tt by rescaling and bringing in an additional numerical 
factor multiplying the angular term, or by also rescaling r and bringing in a numerical 
factor multiplying the dr^ term. 

Using the standard known expressions for the Christoffel symbols, Riemann 
curvature tensor, and Ricci tensor associated with a given metric [71 p 134, 159, 164], 
all of the components of these objects can be calculated for this static, cylindrically 
symmetric metric. The results for this are presented below. 
Nonzero Christoffel Symbols: 

= = $' 
r;^ = A' 

n<p = --e-^^ (3) 

1 



■p<p _ -pip 

rz ZT 



Nonzero Riemann Curvature Tensor Components: 



Static solutions of Einstein's equations with cylindrical symmetry 



4 



R\^^ = -($" + $'2 _ $'A')e2(*-A) (4) 
Nonzero Ricci Tensor Components: 

r 

i?^^ = - $'2 + $'A' + ^A' - - + A'^' 

= r(A' - $' - vl/')e-2A (5) 

R^^ = _(^" + ^'2 _ ^'A' + + I^')e2(*-A) 

r 

Primes correspond to differentiation with respect to r, e.g., $' = 

We would like to solve the Einstein field equations for the vacuum solution, which 
corresponds to Gap = 0. It is easy to show, however, that it is sufficient to calculate 
the solutions for R^p = 0. We begin with the standard definition of the Einstein 
tensor. Gap = Rap — \Rgap- From this we can calculate the trace of the Einstein tensor 
G := G^^ = R^^ — \Rg^^ = R — 2R = —R and thus obtain the following relation between 
the Ricci and Einstein tensors: Rap = Gap — \Ggap- Thus we see that if Rap = then 
Gap = 0, and conversely, if Gap = then Rap = 0. Thus the solutions to Rap = are 
also the solutions to the vacuum Einstein field equations. Gap = 0. 

By equating the nontrivial components of the Ricci tensor with zero, we obtain a 
set of four ordinary differential equations for $, A, and We further note that the 
exponential function is never equal to zero, so the differential equations reduce to 

($" + $'2 _ $'A' + i<l>' + ^'$') = 0, (6) 

r 

(A' - $' - = 0, (8) 

r 

We see that ([8]) can be solved for A' in terms of the other two unknown functions, which 
can then be substituted into ([6]), ([7]), and ([9]) to eliminate A'. Thus this system can be 
reduced to 

A' = $' + (10) 

$" + i$' = 0, (11) 
r 

^" + = 0, (12) 
r 

$'^' + ^$' + ^^' = 0. (13) 
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Now (ITUI) , (ITT]) , and (IT^ are linear second-order equations easily solved by separation of 
variables, yielding $ = ln(r"i) + ln(a2), ^ = ln(r''i) + ln(62), and A = ln(r"i+''i) + ln(c). 
Also, ( fT3|l provides the additional constraint that aibi + ai + 61 = 0. Thus the static, 
cylindrically symmetric metric is 

ds' = -ay^dt' + c'r'(-^+'^)dr^ + rM^^ + b^^dz^ (14) 

with < < 0*. The multiplicative constants 02 and 62 can easily be absorbed by 
rescaling t and z, resulting in 

ds^ = -r^-^dt^ + c^r^i-^+'')dr^ + rM^^ + r^^M^^; (15) 

after each change of variables in what follows, we shall carry out this procedure again 
without comment. Here we have shown that the coefficients must be powers of r as in 
( IT5I) . with aibi + ai + 61 = 0. Since we no longer have to worry about the constants 02 
and 62? "we now drop the subscripts on ai and 61, and simply write the constraint as 

ab + a + b = 0. (16) 

The constant c can also be absorbed by rescaling r, which affects the d0^ term by 
bringing out another constant in front, resulting in 

ds^ = -r^-dt^ + r^(-+'')dr^ + K^'^dcj)^ + r^Mz^. (17) 

This leads to two natural conventions for the dcj)^ term. One can now rescale so that the 
constant K"^ is absorbed, thus redefining the range 0* of 0. One could instead rescale 
to fix its range to be from to 2-7?, in which case the constant remains, multiplying either 
d0^ as in (fTTl) or dr^ as in (fTSl) . In the work that follows, we use the first convention, 

ds^ = -r^-dt^ + r2(«+^)dr2 + rM(j)^ + r^'dz^, (18) 
where the arbitrary constant is hidden in the periodicity, 0*. 



3. Special cases 

We now examine in greater detail the relationship between a and b, which is illustrated 
in figure [H 

We note the existence of several special points on these graphs and examine their 
significance in various different forms of writing the cylindrical metric. One such point 
is a = 6 = 0, which reduces the metric of ( fTSl) to 

ds^ = -dt^ + dr^ + rM(f)^ + dz\ (19) 

This describes a cone; that is, flat space missing a wedge of deflcit angle A0 = 27r — 0*. 
If 0* > 27r, a wedge is added. Ordinary Minkowski space arises as the very special case 
0^, = 2tx. It is also useful to note the symmetry between the t and z coordinates, along 
with the Lorentz symmetry under boosts in the z direction. 

Another point of interest is 6 = —1, in which case a — )■ 00. The signiflcance of this 
(apparently singular) case can be better demonstrated if we rescale r, t, and z of (fTTll 
(with the rescaling for r given explicitly in Section 4) to write the metric in the form 

ds^ = -r-^'dt^ + r2(i+^)d02 + AV2^(i+^)(dr2 + d^^) (20) 
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Figure 1. (a) Plot of the relationship between a and h. (b) Plot of the relationship 
between a and a + 5. 



with A = c{l + b). 

If we now treat A as the arbitrary constant instead of 6, the metric remains 
nonsingular when 6 = — 1 in the other terms. After rescahng r with f = A^^^^'^^^^^^^^^'^ 
to absorb A, one gets 

^^2 ^ _--2fe^^_^^[2(l+fe)]/[-b(l+6)-l]-2(l+6)^02_^-2fe(l+b)^-2_^ -26(1+6)^-2^ (^21) 

We must now rescale as well in order to absorb the final constant in front of the d0^ 
term; this changes the range 0^,. The metric becomes 

= -f-^M? + f2(l+'')d02 ^ ^26(1+6)^-2 ^ ^26(1+6)^-2^ (22) 

and at the point where 6 = — 1 it reduces to 

ds^ = -f^dP + d02 + df2 + df. (23) 

Under the transformation T = fsinht and R = fcosht we see that (l23|) is locally 
equivalent to fiat space, 

ds^ = -dT^ + d02 + dR^ + df (24) 

with a periodic coordinate. 

We also note that the general relationship in ( IT6|) is symmetric when a and 6 are 
switched, corresponding to switching z and t. This observation suggests that the case 
a = —1, 6 — 7- oo is parallel to the foregoing one. To see its physical significance, we can 
write the metric of ( 1T5|) in the form 

ds' = L-2V2-^+2-(-dt2 + dr^) + L2(^+^)r^+2a^^2 ^ ^-2(l+6)^-2a^^2 (25) 

with L = [(1 + b)/K]^^^''^ ^ . At the point where a = —1 and 6 — )• oo, this reduces to 

ds^ = -dt^ + dr^ + d^^ + rW. (26) 
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Under the transformation Z = r sin z and R = r cos z we get 

ds^ = -dt^ + d02 + di?2 + dZ\ (27) 

and this once again looks hke flat space locally, but with still a periodic coordinate. 
We note that these locally fiat solutions are not included in the general solution found 
in Section 2 because there we fixed the coefficient of d0^ to be r^, whereas in and 
( 1271) that coefficient is a constant. 

4. Transforming between metric conventions 

The equivalent forms ( IT5|) . (!20|) . and ( 125|) all appear in [8j. The historical rationale for 
the last two will be explained in Section 6. 

We now examine how transforming from ( fTSl) to the other metric forms affects the 
radial coordinate r. To go from f|T5|) to fl20|) . we must use f = [c/{a + l)]r'^+^. We 
see that in this case, the exponent of r is negative whenever a < —1 or, equivalently, 
6 < — 1. Under this condition, r = in our gauge choice corresponds to f = oo. To 
go from our form of the metric to that of ( l25l) . we require f = [c/{b + l)]r^+^. Once 
again, the exponent of r is negative whenever b < —1 {a < —1), and in that case r = 
corresponds to f = cxd and vice versa. 

Another natural choice for writing the metric, which we shall call "arc-length 
gauge," is 

ds^ = -Adt^ + Bd(f)^ + dr^ + Cdz\ (28) 

where A, B, and C are once again functions of r only. In this case, to transform from 
our gauge to fl28l) we require f = [c/{a + b + i"jyo.+b+i^ Under this transformation, the 
metric becomes 

= - (Df) + + (Df) [2/('^+''+l)] + ^Df) [2V(a+6+l)] d^2 (39) 

with D = {a + b + l)/c. The exponent of r in our definition of f is negative whenever 
a + b < —1, which occurs whenever b < —1 {a < —1). Thus in all three alternate 
metric forms discussed here, r = in our gauge corresponds to f = oo in the new 
gauge whenever a < —1 and 6 < — 1 (hence a + b < —4 from figure [1] (b)). The other 
possibilities have a + b > and r and f running in the same direction. 

It is also interesting to calculate W = R""^^^ Raj3iJLv for the vacuum solution, because 
this is the simplest nonzero curvature invariant (since i? = from G^/j = 0). We get 
the result that 

R^^^'Ra^f.u = 4Cr-^("+*+i) (30) 

where C = 3a^ + 3b^ + 3a^6^ + 2ab + 2a% + 2a6^. We see that the exponent is negative 
whenever a + b > —1, and in that case, — )■ as r — )• oo. In the other case, W ^ oo 
as r — )■ oo but — )■ as f — oo, and in that case there is some ambiguity as to which 
of these limits is the "outside" and which is the "inside" . 
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5. Solutions of the Einstein equations with sources 

We would like to find some cylindrical space-times that are not singular along the 
central axis. This requires solving the Einstein equations in cases where T has nonzero 
components. In order to proceed with this, we first require a few more basic quantities 
and tensors encountered in general relativity. We present here the results for the Ricci 
scalar, i?, the Einstein tensor, G^^, and the stress-energy tensor, T^j,, for the cylindrically 
symmetric metric given in ([2]). 
Ricci Scalar: 



2 2 2 

f A' _ f $' _ f (31) 

/y /y* /y* 



R = e-^^{-2^" - 2^'^ + 2<I>'A - 2^" - 2^'^ + 2^'A' - 2^'$' + 
2 2 2 

/y /y 

Nonzero Einstein Tensor Components: 

Gu = e2(*-^)(-^" - ^'2 _^ ^/A' + 1a' - ^^') 

Grr = + + 

G^^ = r\-^^{^" + - $'A' + ^" + - ^'A' + ^'$') (32) 

G,, = e2(*-A)($" + $'2 - $'A' - ^A' + ^<l>') 

The components of the stress-energy tensor are defined by setting, for example, T'^r = Pr 
and lowering an index to get Tj-r = Pr^^^- The other pressure components are defined 
similarly, and T\ = —p. 
Nonzero Stress Tensor Components: 



Tu = pe 




rp 2A 




T^<l, = P<i>r'^ 


(33) 


Tzz = PzQ 





From the Einstein field equations in natural units, G^j, = SttT^;^, and the conservation 
laws, T"^;;3 = 0, we get the following five differential equations: 

= ^ + Pr{^' + ^'+^-) + p<f ' - Pz^' - ^p^, (34) 

47r(p + pr+p^ + Pz)e^^ = + - $'A' + + -<^', (35) 



r 



4:7l{p + Pr -P4>-Pz)e 



2A 



_ $'2 + $'A' + -A' - - + A'^', (36) 

r 

47r(p - +p^- p,)e2^ = 1(A' - $' - M/'), (37) 

r 

47r(p -Pr-P4>+ Pz)e^^ = - + ^'A' - - (38) 
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We can simplify these by summing and and subtracting (155]) . This yields 

47r(p + 3pr +p^- p,)e^^ = 2^'$' + -<l>' + + -A'. (39) 

We now add and subtract equation ( 1371) from ( l39l) . resulting in 

47r(2p + + 2p0 - 2p^)e^^ = 2^'$' + ^A' (40) 

47r(4p,)e2^ = 2^'$' + -($' + ^')- (41) 

r 

We now have a system of differential equations where equations (1M|) . ( 15^ . (157|) . and 
( 138|) can be solved for p^, Pz, ^, and A (given p and an equation of state relating p 
and the various pressures), and (HT]) . which contains only lower-order derivatives of the 
unknown functions, provides an additional constraint. The system of all five equations 
is second-order in $ and \l/ and first-order in A and pr- 

Differentiating equation (HTj) with respect to r and using equations (1M|) . ( 15B]) . ( 155]) . 
and (j^Ol) to substitute for pj.. A', and \E'" yields an expression which reduces to = 0; 
thus equation (HTl) must hold for all r if it holds at any r. 

5.1. Solutions with p = —pz, Pr = Pel) = 

Solving these equations for arbitrary p, pr, Pcf,, and Pz is rather difficult, so a simpler 
case one can look at is when p = —pz and the other pressure components are zero. In 
this case, the differential equations reduce to 

= p($' + ^'), (42) 

= $" + $'2 - $'A' + ^'<l>' + (43) 

r 

47r(2p)e2^ = -(A'-<l>'-^')> (44) 
r 

= - ^'^ + ^'A' - - (45) 

r 

= 2$'^' + -($' + ^')- (46) 
r 



From equation ( H2l) we can see that $' -|- \E'' = 0, allowing us to solve equation ( l44l) 
for A' = Snpre^^, which can easily be solved using integration by parts to get A for a 
given p. From (j46]) and the fact that $' \1'' = 0, we also see that = 0. Thus 
we can conclude that $' = \E'' = 0, yielding $ = ai and \E' = 02 (where Oi and 02 are 
constants). The metric of the solution can be written as 

ds^ = -dt^ + e^'^dr^ + r^d^^ + dz^ . (47) 

This solution (with p = 1/ {Stttq) where tq is a constant) is widely known; it is usually 
called the "cosmic string solution" or "Gott's solution" |11[5]. Using this value of p, our 
solutions yield 

ds^ = -dt^ + [rl/{rl - r^)]dr^ + rM^^ + dz\ (48) 
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which agrees with |3] after making the substitution r = sin(^) and then rescahng 
coordinates as necessary. This metric is Lorentz-invariant under boosts in the {z,t) 
plane, and thus p and pz are independent of frame. If we did not have the condition 
that p = —pz, the solution would not be Lorentz invariant in this way and p and Pz 
would not be frame-independent. Indeed, generically one would expect the density of 
the matter in a string source to be affected by Lorentz contraction when one moves out of 
the rest frame. But a Gott string, like cosmological dark energy (where all components 
of p equal — p), has no preferred rest frame. 

5.2. Numerical solutions 

We now present some numerical solutions (calculated with Mathematica) for the case 
when p is constant out to a radius R and zero outside of this radius, and the pressure 
is isotropic in all directions, pr = p^j, = Pz = P- This is analogous to isotropic pressure 
in the spherically symmetric case. Since our differential equations involve factors of 
1/r, they present problems when trying to solve the system numerically starting from 
r = 0. In order to deal with this, we first make power series expansions of p, 
and A around r = 0. We keep terms up to order r in the p and A expansions (since 
our differential equations involve first-order derivatives of these functions) and keep 
terms up to order in the $ and \l/ expansions (since the differential equations involve 
second-order derivatives of these functions), resulting in 



P = Po + Pir, (49) 

A = Ao + Air, (50) 

$ = $0 + $ir + $2r^, (51) 

^ = ^0 + ^^r + ^'2^-^ (52) 



We also take the initial conditions \E' = $ = A = 0atr = 0so that the corresponding 
metric coefficients are equal to 1 at that point, and choose = $' = to get smooth 
solutions at the axis. Equations ( H9l) -( l52l) should satisfy our differential equations 
near r = 0, so we substitute them into flM|) . f l35|) . f l37|) . fl38|) . and fHT]) (taking 
Pr = P<i> = Pz = p), and obtain the relationships pi = 0, Ai = 0, $2 = 7r(p + 3po), 
and \l/2 = — 7r(p ~ Po). After choosing values for p and po = p{0), we can determine the 
values of p, $, and A at some small r away from 0; we take r = 0.01. We use these 
as our initial conditions for the numerical calculations and obtain solutions for various 
values of p and po! two examples of such solutions are provided. For the case of p = 1, 
Po = 0.1, the results are given in figure El When p = 10, po = 1, the results are given in 
figure 131 In both cases, the numerical solutions we obtain are not conformally flat (i.e., 
there are nonzero components of the Weyl tensor). 

5.3. Connecting to exterior solution 

After finding interior solutions numerically for Pr = P(f> = Pz = P, we can then connect 
them to the exterior vacuum solution found in Section 2. We take R to be the point 



Static solutions of Einstein's equations with cylindrical symmetry 



11 




Figure 3. For p = 10 and po = 1 (interior): (a) Plot of A(r) (thick), 'i'(r) (normal), 
and <i>(r) (dashed), (b) Plot of p(r). 



where p{r) = 0, and use the values of $(-R), $'(-R), "^{R), "^'{R), and A(i?) from the 
numerical solutions as conditions to determine the unknown coefficients (ai, 02, &i, 62, 
and c) of $, \E', and A from the vacuum case. We must match our interior solution 
with the most general vacuum solution from Section 2 which includes all the arbitrary 
constants, since we chose our initial conditions so that the interior metric coefficients 
are all 1 at r = 0. Because of this, we are not free to scale away 02, 62, and c, and we 
must keep them in the metric in order to match our two sets of solutions. 

To perform the matching, we impose the conditions that $(-R), ^'{R), "^{R), \E''(i?), 
and A(i?) must be continuous at the boundary where p = 0. These five continuity 
conditions provide us with the information needed to find values for the five arbitrary 
constants from the vacuum solutions. After matching the solutions in this manner, A'{R) 
is not necessarily continuous at the boundary. To explain this, it is helpful to look at the 
metric in the arc-length gauge, given in (|28l) . Obviously one wants A, B, and C to be 
continuous so that the metric is continuous. Since r has the same geometrical meaning 
on both sides of the surface, a nonsingular metric should also have continuous A', B', 
and C (where primes indicate differentiation with respect to the radial coordinate in 
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the arc-length gauge). Technically, this condition is called "continuity of the second 
fundamental form" [9]. Also, in this case, the coefficient of dr^ is unity, so it and its 
derivative are automatically continuous. Now letting f be the radial coordinate in arc- 
length gauge, ( 128|1 . and letting r be the radial coordinate in tangential gauge, ([2]), we 
can see that ^ = e^^. Hence, after imposing the conditions that A, $, and \1/ must 
be continuous, continuity of the second fundamental form in arc-length gauge is merely 
equivalent to continuity of ^ and ^ in tangential gauge, and the derivative of the 
coefficent of d(f)'^ in tangential gauge is automatically continuous as well. Thus there are 
no conditions which require that ^ also be continuous at the boundary. 

The results for the two sample cases given above are presented in the figures. When 
p = 1, Po = 0.1, we get that R = 0.1486, and the resulting exterior solutions are plotted 
in figure HI For the case of p = 10, po = 1, we get R = 0.04791, and the exterior 
solutions are plotted in figure O Numerical constants for these solutions are given in 
table [H Also, after rescaling coordinates appropriately in order to put the metric in 
the form of equation fllSI) . the range of changes as well; it had to be 27r for the inner 
solution to guarantee smoothness at the origin, and thus the outer solution initially 
has range 27i when matched with the inner solution. The new value of 0* is given by 
0^, = (27r)c^^/*-"^+^^+"'^-'. The values of 0* for the two solutions described above are also 
given in the table. 

Bicak et al. have done a more general study of static perfect-fiuid cylinders, 
including numerical work on incompressible cylinders and their external vacuum 
solutions. We have compared our numerical results with theirs and found that they 
agree fairly well. We made comparisons with their parameters m and C, which are 
related to our notation by m = —bi and C = {2n/(f)^:)(ai + l)[^i/(''i+^i+i)]. The results 
are presented in table [2] and table [31 

In the spherically symmetric case, the Buchdahl theorem [71 p 269] requires that 
R > |M for any stellar model, where M = ^npR^. This implies that R^p < ^, 
where R^p is a dimensionless quantity. Although this theorem does not apply to the 
cylindrically symmetric case, it is interesting to note that the inequality does hold in 
the examples studied above. 
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Figure 4. For p — 1 and po = 0.1 (exterior): Plot of A(r) (thick), 5'(r) (normal), and 
$(r) (dashed). 
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Figure 5. For p = 10 and po ^ 1 (exterior): Plot of A(r) (tlrick), ^'(r) (normal), and 
<l)(r) (dashed). 

Table 1. Numerical constants for exterior solution when p — I and po =0.1 (second 
column), and p = 10 and pq — 1 (third column). 



R 


0.1486 


0.047 91 


R^p 


0.022 08 


0.022 96 


ai 


0.2052 


0.2136 


h 


-0.1703 


-0.1761 




1.627 


2.114 


b2 


0.6722 


0.5435 


c 


1.279 


1.342 


(t>* 


4.954 


4.732 



Table 2. Comparison of our numerical results with those of Bicak et al. fTO], for p = 1 
and po = 0.1. 

Bicak et al. our results 
TO 0.1696 0.1703 
C 1.05901 1.05901 



Table 3. Comparison of our numerical results with those of Bicak et al. [TOj, for 

p = 10 and po = 1. 

Bicak et al. our results 
TO 0.1696 0.1761 
C 1.10063 1.10179 



6. History 

The static, cylindrically symmetric vacuum solutions were found in the early 20th 
century by Weyl [2] and Levi-Civita [3]. They were interested in the more general 
problem of static geometries that are merely axially symmetric, that is, depend on z as 
well as r. It was therefore convenient to use coordinates that treat r and z on the same 
footing, so the Weyl-Levi-Civita solution was found in the form (|20|) . 
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Similarly, Marder [H], building on work of Rosen [TT], studied cylindrically 
symmetric gravitational waves and hence treated r and t alike, obtaining the static 
solution in the form ( l25|l (four decades after [2113]). Marder's paper, which also displays 
versions (fTSjl and ( 120|) . was very helpful in the analysis of vacuum solutions in the present 
paper. 

The arc- length gauge (128]) . leading to solution fl29|l . was probably first used by 
Evans [12]; it appears also in much of the cosmic-string literature, such as [6], although 
the Weyl-Levi-Civita convention is also popular there. Of course, there are many other 
possible gauge choices; for example, Bicak et al. [10] use the convention (in our notation 
(CD) $ = -A. 

During the 1980s, cylindrical solutions received broad attention when the (locally 
flat) cone solutions were studied as representing the space-time outside a one- 
dimensional concentration of matter or gauge-field energy, a "cosmic string". Classic 
papers by Taub [13] and Vilenkin [T3] considered infinitely thin sources; nonsingular 
solutions with sources of finite radius were constructed by Gott [1] and Hiscock [5] ; and 
both were studied by many more physicists. With few exceptions (e.g., [I5l|16]), the 
possibility that the exterior of a cosmic string might be one of the nonflat Weyl-Levi- 
Civita solutions was not widely recognized, except when the theory was generalized 
to include a cosmological constant [171[T8] (more recently studied in [19]). Physical 
enthusiasm for cosmic strings as realistic cosmological objects has diminished in recent 
years under the pressure of new observational data. 

Nevertheless, research on static cylindrical solutions with interior sources and 
general Weyl-Levi-Civita exteriors has continued and intensified in recent years 
[9l [T0l[T2l[T9]434] . We are not able to provide a complete review of the literature here. 

7. Conclusion 

Cylindrical symmetry in general relativity turns out to be similar to spherical symmetry 
in many ways but quite different in others. The only static vacuum spherical solutions 
are the Schwarzschild metrics parametrized by mass, but, contrary to popular belief, 
there are static vacuum cylindrical solutions other than the cone, or cosmic string, spaces 
parametrized by defect angle. A cone solution is Lorentz-invariant along the cylinder 
axis and hence cannot arise from a matter source (at least in the small-radius limit) 
unless the latter has a very unusual equation of state. We have presented numerical 
solutions with more conventional interior sources and more general exterior geometry. 

The construction of the various solutions presented here illustrates several 
instructive points. The choice of gauge (coordinate system) is always a major issue 
in relativity; the same space-time can look quite different in different gauges, and how 
(and whether) to choose a standard gauge or "normal form" for a given problem is not 
always obvious. For our problem there are several natural ways to fix the gauge, and we 
have taken pains to describe them all and how they are related. Even after a definition 
of radial coordinate has been selected, further steps to a normal form can be taken by 
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linear rescaling of the coordinates. But, as our final calculation shows, sometimes this 
progress must be undone to match solutions for different regions properly. 

The structure of the Einstein equation system is nontrivial. There is one more 
equation than one might naively expect. The extra equation serves as a constraint 
on the data. For the cylindrical vacuum solutions this constraint is a simple algebraic 
relation among the parameters, but for our interior solutions, extra work was needed to 
verify that the constraint equation is consistent. The differential orders of the remaining 
equations need to be considered carefully in order to choose the correct sort of initial 
data on the axis and to match the interior solution properly to an exterior vacuum 
solution. 

Finally, we observed some surprising ambiguities of interpretation. On the left 
branch of the vacuum solutions (see either half of figure [1]) an increasing radial metric 
component corresponds to decreasing behavior of some of the other metric components; 
therefore, what is naturally considered the radial coordinate in our gauge is naturally 
considered to be the reciprocal of such a coordinate in other gauges. (That is, our 
axis is previous authors' infinity.) Also, two special solutions, f|T9|) and f l27|) . are not 
really "cylindrical" spaces, but rather representations of flat space-times in nonstandard 
coordinate systems. 
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